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Calibration of three scintillators (EJ232Q, BC422Q and EJ410) in a time-of-flight (TOE) arrangement using 
a laser drive-neutron source is presented. The three plastic scintillator detectors were calibrated with gamma 
insensitive bubble detector spectrometers, which were absolutely calibrated over a wide range of neutron 
energies ranging from sub MeV to 20 MeV. A typical set of data obtained simultaneously by the detectors 
are shown, measuring the neutron spectrum emitted from a petawatt laser irradiated thin foil. 


I. INTRODUCTION 

In recent years, there has been a growing interest in 
laser-driven sources for their high brightness, direction¬ 
ality and compactness, which have shed a new light on 
the production of ultra-short (sub-ns) neutron pulses that 
can be used in many potential applications in science^, 
industry^, security^ and healthcare^’^. Since laser-driven 
neutron sources are primarily based on beam-fusion or 
spallation reactions initiated by fast ions^, the measure¬ 
ment of absolute neutron spectrum is crucial not only for 
the development and optimisation of the neutron sources, 
but also for the study of the parent ions involved in the 
neutron generation process. 

Among commonly used neutron diagnostics, those 
based on time-of-flight (TOE) arrangement are proba¬ 
bly the most widely used for characterising the neutron 
spectra. Depending on the distance between the source 
and the detector, a TOE spectrometer can in princi¬ 
ple resolve neutron energies ranging from sub-MeV to 
hundreds of MeV. However, the most unique feature of 
TOE scintillator detector lies in the capability to pro¬ 
vide an absolute spectral measurements of fast neutrons, 
which are otherwise first moderated to lower energies in 
other types of neutron spectrometers resulting in a poor 
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spectral information^. Most common TOE detectors are 
made of scintillators coupled to photo-multiplier tubes 
(PMT), where the specifications of the scintillators and 
PMTs depend on the type and range of measurements 
sought. Plastic scintillators are the most common type 
of detectors employed for neutron detection due to their 
high dynamic range and ease of handling. 

In order to obtain the complete spectral information 
from a TOE diagnostic, one would require an abso¬ 
lute calibration of the detector over a range of neu¬ 
tron energies. Traditionally, the scintillators are cali¬ 
brated with continuous sources, such as ^^^Cf emitting 
neutrons via spontaneous fission reactions^’^, by using 
pulse-height discrimination technique in multi-channel 
analyserHowever, the radioactive isotopes do not pro¬ 
vide an accurate measurement of detector response due 
to their broad energy spectrum (for example, 0.1 - 10 
MeV for 

Since the scintillators fluoresce by the knock-on pro¬ 
tons produced by neutron elastic scattering, their re¬ 
sponse to different energy neutrons exhibit a highly non¬ 
linear behaviour. In this paper we report on an alterna¬ 
tive way of calibrating the scintillator detectors by using 
an ultra-short burst of neutrons, which enables an en¬ 
ergy resolved calibration by deploying the detector in a 
TOE arrangement. High power laser driven pulsed neu¬ 
tron sources emit high flux of neutrons within a short 
burst duration (sub-ns). Therefore using a detector at a 
suitably large distance (several meters) from the source 
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FIG. 1. Schematic of the experimental setup. The neutron beam generated from laser irradiated target or catcher propagates 
through various diagnostics before reaching the scintillator detectors (Figure not drawn to scale). The scintillator detectors 
were placed face-on to each other inside a lead and plastic housing. The collimators, made of 5 cm thick plastic, were used 
before the detectors in order to reduce stray neutrons hitting the detector. 


temporally disperses the arrival of different energy neu¬ 
trons on the detector. This provides an ideal situation 
for energy resolved calibration of the detector using a 
broadband source. This is unlike the pulse-height dis¬ 
crimination technique, where the neutron energy is iden¬ 
tified from the height of the pulses generated by neutron 
induced scintillations, while assuming that each pulse is 
produced by a single hit of neutron. The laser driven 
source used in our case delivered high flux of neutrons 
of the order of 10 ^ n/sr in a single shot, which was able 
to produce a strong signal in the TOF detector for cross¬ 
calibration with the absolutely calibrated bubble detector 
spectrometers (BDS) fielded in the same shot. 

The calibration technique presented here combines a 
series of experimental measurements and Monte-Carlo 
simulations, which was employed to calibrate three de¬ 
tectors made of different plastic scintillators. As a typi¬ 
cal detector consists of many parts, such as scintillator, 
light guide and photo-multiplier tube (PMT), an abso¬ 
lute calibration requires specific knowledge of the detec¬ 
tor design and characteristics of its different parts, such 
as the detector housing, light guide, PMT etc. However, 
the non-linear response of the detector to different energy 
neutrons comes from the plastic scintillator itself, which 
originates from the difference in recoil proton spectra pro¬ 
duced by different energy neutrons. The relative response 
of the scintillator over a range of neutron energy was ob¬ 
tained by monte-carlo simulations, while benchmarking 
the results with the literature. The neutron spectra from 
the TOF data was obtained by using the energy depen¬ 
dent scintillator response, while taking into account the 
transmission of different energy neutrons through various 
objects along the line of sight to the detector. Finally, the 
detectors were calibrated by comparing the TOF neutron 
spectra with that obtained by the absolutely calibrated 
BDS spectrometers employed in the same shot. 

The subject matter in this paper is organised in four 
sections. Section II illustrates a brief description of the 
experimental setup. Section III describes the basic prin¬ 
ciples of neutron detection by scintillation processes. Sec¬ 
tion IV outlines the numerical work carried out by em¬ 


ploying Monte-Carlo simulations in order to obtain the 
relative response of our scintillator to different energy 
of fast neutrons. Section V reports the final calibration 
obtained for the three scintillator detectors and typical 
neutron spectra obtained by different detectors. 

II. EXPERIMENTAL SETUP 

Calibration of the scintillator detectors was carried out 
using a laser driven neutron beam generated at VUL¬ 
CAN Petawatt target area of Rutherford Appleton Lab¬ 
oratory, STFC, UK^i. The 750fs FWHM laser pulse of 
Vulcan, with energy of 600J, was focused on the target 
by a f/3 off-axis parabola delivering intensity in the range 
of 1 — 5 X lO^^Wcm~^ . Energetic ions were produced 
by irradiating the laser on deuterated plastic foils of vari¬ 
ous thickness, which were employed to generate neutrons 
through beam fusion reaction by interacting with either 
the bulk target ions or secondary C 2 D 4 catcher target 
placed at a small distance (^cm) from the laser irradi¬ 
ated target. 

A suite of diagnostics, such as scintillator TOF, BDS, 
activation samples^^’^^ and CR39 stacks^^~^^ were em¬ 
ployed around the target in order to characterise the neu¬ 
tron beam parameters. A schematic of a typical exper¬ 
imental setup used for cross-calibrating the scintillator 
detectors with the bubble detectors is shown in Fig. 1 . 
The bubble detectors, manufactured and calibrated by 
BTI Bubble technology industries'^, were placed inside 
the interaction chamber at various distances (30-50 cm) 
from the neutron source in order to ensure a sufficient 
number (at least ^ 50 bubbles per energy bin) of bubble 
formation. The scintillator detectors were placed outside 
the interaction chamber at a significantly larger distance 
from the neutron source to achieve high energy resolu¬ 
tion by the TOF technique. The scintillator detectors 
were shielded appropriately by lead and plastic blocks for 
the purpose of reducing noise level mainly originated by 
gamma rays and thermal neutrons hitting the scintilla¬ 
tors. Two 5cm thick, one-meter long plastic collimators. 
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with 10 cm X 10 cm internal aperture, were used along 
the line of sight of the scintillators in order to reduce 
stray neutrons hitting the detector. 

The three plastic scintillators discussed in this paper 
are EJ232Q^^ (0.5% quenching, 150 mm diameter, 25 
mm thick), BC422Q^^ (1% quenching, 180 mm diam¬ 
eter, 10 mm thick) and EJ410^^ (90 mm diameter, 16 
mm thick). Calibration of EJ232Q and EJ410 scintilla¬ 
tors was obtained in the same set of shots while using 
the two detectors inside the detector housing as shown 
in the Eig. 1. However, the BC422Q scintillator was cal¬ 
ibrated with the BBS in a separate set of shots where 
the BC422Q detector was kept alone inside the detector 
housing. Where EJ232Q and BC422Q have fast rise and 
decay times (sub ns and a few ns respectively), EJ410 has 
significantly slower decay of 200 ns. Each scintillator was 
coupled with a fast PMT tubes for converting the light 
output from the scintillator into electrical signal as well 
as to amplify the signal for detection. The gain of PMT 
was controlled by varying the biasing voltage from shot 
to shot in order to avoid signal saturation. The detectors 
were attached to fast (6GHz) oscilloscopes for recording 
the output from the PMTs. 


III. BASICS OF SCINTILLATION PROCESS 


Since neutrons do not interact via the Coulomb force, 
they don’t produce excitations in the scintillation ma¬ 
terial directly. Therefore neutron detection by scintilla¬ 
tors relies on secondary ionising particle generation in 
the scintillator material via elastic, inelastic or nuclear 
reactions. East neutrons (generally ^0.5 MeV) can ef¬ 
ficiently produce recoil protons via elastic scattering in 
hydrogen-rich materials, whereas slow neutrons rely on 
nuclear reactions such as the (n, 7 ) or (n,(a) reactions to 
produce ionisation. Materials such as Lil(Eu) or glass sil¬ 
icates doped with ^Li and ^^B are therefore particularly 
well-suited for the detection of slow (thermal) neutrons. 

The response of plastic scintillators to fast neutrons de¬ 
pends on several factors. The dominant mechanism for 
production of secondary ionising radiation for incident 
neutron energy up to a few 10s of MeV is the proton 
recoil based on elastic n-p scattering. As neutrons of 
a given energy pass through the scintillator, they pro¬ 
duces a spectrum of recoiled protons by direct elastic 
scattering with hydrogen atoms present inside the mate¬ 
rial. Energy of the recoiled protons varies from zero to a 
maximum value, which is equal to the incident neutron 
energy. Each recoiled proton then produces scintillation 
light, Lp{E), while travelling in the scintillator over a fi¬ 
nite length depending on its kinetic energy. A part of the 
energy deposited by the protons is used for the scintilla¬ 
tion process and the Lp{E) can be calculated with semi- 
empirical formulae described by Chou^^ and Wright^^. 


Where Chou’s formula, given as 


dLp{E) ^ ^dE 
dx dx 


1 + KB^ + C 

dx 


\ dx 


( 1 ) 


which is a generalisation of the Birk’s formula^^, 
Wright formula. 




dx. 


( 2 ) 


on the other hand uses a simpler model while attribut¬ 
ing for quenching effects in scintillators. In the above 
formulae, *5, KB^ C, A and a are the fitting parameters 
to the experimental data, R is range of the proton inside 
the scintillator and dE/dx is the proton energy loss in 
the scintillator material per unit length. 


IV. RELATIVE RESPONSE OF EJ232Q TO FAST 
NEUTRONS 

Although Chou’s and Wright’s semi-empirical formulae 
provide fairly accurate estimation of the light output, the 
statistical nature of particle scattering and energy loss 
mechanisms demands a thorough monte-carlo approach 
in order to characterise the response of plastic scintilla¬ 
tors to fast neutrons. Most of the commercially available 
plastic scintillators are experimentally calibrated for dif¬ 
ferent types of ionising radiation, such as electrons, pro¬ 
tons, alpha particles etc., as provided by the manufac¬ 
turers. However the calibration of plastic scintillators for 
fast neutrons is often unknown, as is the case for EJ232Q 
and BC422Q used in our experiment. Therefore we used 
a simple model based on a series of systematic monte- 
carlo simulations, as described below, in order to obtain 
a relative response in terms of the light output by the 
whole spectrum of recoiled protons. 

At first the energy spectrum of recoiled protons pro¬ 
duced by a given energy of neutrons in EJ232Q plas¬ 
tic (which is the same material used in BC422Q) was 
obtained by employing ELUKA^^ simulations. The sto¬ 
ichiometric information about the EJ232Q plastic was 
taken from the manufacturer documentation. A series of 
ELUKA simulations was done by irradiating a small cube 
of the scintillator material (lOqm x lOqm x 10pm) with 
different neutron energies. Eor each neutron energy the 
recoiled protons generated inside the cube were tallied. 
As an example, Eig. 2a shows the recoiled proton spectra 
produced by neutrons of three different energies. 

The next step of calculation involves finding the stop¬ 
ping range of protons in the EJ232Q material, which was 
simulated using SRIM^"^ monte-carlo simulation. The 
stopping range was calculated by obtaining dE/dx as 
a function of incident energy of protons while passing 
through a very thin (1 pm) layer of EJ232Q, as shown 
in Eig. 2b. Using the stopping range of protons, Lp{E) 
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FIG. 2. a) Spectra of recoiled protons generated inside a 10pm x 10pm x 10pm cube of EJ232Q scintillator for different incident 
neutron energies, b) shows stopping range of protons inside EJ232Q scintillation material obtained by SRIM simulation, c) 
shows scintillation light output for different proton energies calculated for EJ232Q scintillator by using Chou’s and Wright’s 
formulas, as well as for a few different scintillators obtained from literature^^. (d) shows the relative, weighted average light 
output from the EJ232Q obtained by using Eq.3, where the recoiled proton spectrum is shown in (a) and Lp{E) is shown in 
(e). 


for the EJ232Q scintillator material was calculated for 
different proton energies by using Chou’s and Wright’s 
semi-empirical formulae, as shown in Fig. 2c. As one can 
see, the Lp{E) from both formulae closely agree with each 
other and with the data provided by the manufacturer. 
It also closely matches with the data found in literature 
for similar plastic scintillators^^. 

Finally, the weighted average scintillation 

light output for a given neutron energy En is obtained 
by integrating the light output from all the recoiled pro¬ 
tons generated by the fast neutrons inside the scintillator. 
Mathematically one can write 


L 


n 



dNp 

He 


Lp{E) dE 

En 


( 3 ) 


where dNpjdE is the recoil proton spectrum shown in 
Fig. 2a and Lp{E) is the light output by different energy 
protons in the scintillator, shown in Fig. 2c. By following 
the steps mentioned above, the weighted average light 
output from the scintillator was obtained for different 
incident neutron energies as shown in Fig. 2d. 

Since the BC422Q is the same scintillator material as 
EJ232Q, produced by a different company, we used the 
same detector response for EJ232Q and BC422Q in the 
data analysis described below. However, for the EJ410 
scintillator we used the relative response of the scintilla¬ 
tor to different neutron energy provided by the manufac¬ 


turer^^. 



FIG. 3. (a) Typical raw data obtained by the two differ¬ 

ent scintillators, EJ232Q and EJ410, in the same shot. The 
sharper peak in the signal near to zero time corresponds to 
prompt gamma irradiation, while the following peak corre¬ 
sponds to neutron signal, (b) shows deconvolved neutron 
spectra from the raw signal shown in (a), after accounting for 
PMT gain and scintillator detection response. 
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V. CALIBRATION OF THE SCINTILLATOR 
DETECTORS 

The TOF arrangement allows us to distinguish the sig¬ 
nal produced by the neutrons from the gamma-rays, pro¬ 
duced during the laser interaction, by their difference in 
arrival time at the detector. The prompt gamma flash in 
the raw data at early time, as can be seen in the Fig. 3a, 
serves as a precise temporal fiducial to identify the time 
of neutron generation at the source. Using this tempo¬ 
ral reference, the time axis of the raw data was con¬ 
verted into neutron energy by taking into account the 
distance between the source and the detector. Due to 
the large distance (^11 meter) between the source and 
the detector employed in the experiment, the signals due 
to gamma radiations and neutrons are well separated in 
time, as can be seen in Fig. 3a, which allowed us to dis¬ 
card unambiguously the gamma signal from the raw data 
before processing further. 

The scintillator detectors used in the experiment were 
coupled with fast photo-multiplier tubes in order to con¬ 
vert the optical photons from the scintillator into electri¬ 
cal signal with a gain commensurate to the applied bias 
voltage. The gain characteristics of the PMTs are cali¬ 
brated precisely by the manufacturers, which were used 
in the data analysis, following the background subtrac¬ 
tion of the raw data and deconvolution of the signal by 
scintillation temporal profile. Assuming the light col¬ 
lection by detector geometry, transmission through light 
guide, quantum efficiency of the PMT etc. are constants 
for a given detector, each point in the processed data now 
corresponds to the light emitted by the scintillator due 
to the neutrons within the respective energy bin. Since 
the energy bins are typically very small (for example, 
^IkeV at 1 MeV and ~10 keV at 10 MeV), the rela¬ 
tive number of neutrons in a given energy bin, dN{En), 
is now obtained by dividing each point of the data by 
the corresponding Ln{En) (given in Fig. 2(d)). Finally, 
the neutron spectrum (dN/dE) was obtained by dividing 
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FIG. 4. Transmission in % of different energy neutrons from 
source to the detector obtained by FLUKA simulations for 
the setup shown in Fig. 1. 


each data point by the width of the respective energy 
bin. Fig. 3b shows neutron spectrum obtained from the 
data collected by the EJ232Q and EJ410 detectors. How¬ 
ever, the spectrum does not yield the absolute number 
of neutrons hitting the detector because of the various 
unknown factors, such as light collection by the detector 
housing, attenuation of scintillation photon by the light 
guide, quantum efficiency of the PMT photocathode etc. 
Therefore, it is important to cross-calibrate the detectors 
against a standard dosimeter. In our case, absolute flux 
calibration was carried out against BDS fielded in the 
same shot along the same direction of observation. 

Bubble detectors (BD) are passive detection devices 
that provide neutron response proportional to the num¬ 
ber of bubbles formed in the detector^^. The response 
of the bubble detectors depends on the kinetic energy 
of the incident neutrons and the composition of the BD 
medium^^. Therefore the BDS consists of six sets of bub¬ 
ble detectors calibrated for different energy ranges of neu¬ 
trons. After exposure to the broadband neutron source, 
the number of bubbles in each of the bubble detector was 
manually counted and the neutron spectrum is obtained 
by deconvolving the number of bubbles in each set of de¬ 
tectors against their calibration curves provided by the 
manufacturer. 

In order to cross-calibrate the scintillator detectors 
against the data obtained in the BDS, it is important 
to assess the transmission and scattering of neutrons by 
the various objects (including the BDS) they encounter 
along the scintillators’ line of sight. Since the scintillators 
are placed at a significantly large distance (^11 meter) 
from the BDS, any small angle scattering of neutrons 
by the intermediate objects would prohibit the scattered 
neutrons from reaching the scintillator detector. On this 
basis the neutron propagation through the system was 
simulated using ELUKA in order to estimate transmis¬ 
sion to the detector as a function of neutron energy, as 
shown in the figure 4. 

In the simulation, the neutron source was defined as a 
pencil beam with no divergence. The simulation was set 
by considering all different objects on the line of sight to 
the scintillator detector, including the chamber wall and 
the detector housing. A series of simulations were carried 
out for different neutron energies at the source and the 
neutron spectra at the corresponding detector was tal¬ 
lied. Due to the heavy shielding (typically 5 cm of lead 
and 5 cm of plastic) used around the detectors (to sup¬ 
press the noise level in the detector by gamma rays and 
thermal neutrons) the scattered neutrons in the vicin¬ 
ity of the detectors can in principle contribute towards 
the signal for that energy in a TOE mode. In order to 
account for the scattered neutrons, two sets of simula¬ 
tions were done for each energy, one with and one with¬ 
out the lead and plastic shielding around the detector. 
The neutron spectra at the detector obtained from these 
two simulations were subtracted in order to estimate the 
spectra of scattered neutrons reaching the detector due 
to the detector housing. The transmission factor is then 
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FIG. 5. (a) Calibrated neutron spectra obtained from EJ232Q and EJ410 detectors compared with that obtained by BDS in 
the same shot, (b) Comparison between neutron flux obtained from different scintillator detectors and the BDS for neutron 
energy in the range 2.5-20 MeV. The graph plots the comparison for three different scintillator detector, EJ232Q, BC422Q and 
EJ410, obtained over a number of shots. The solid lines show the linear fit to the data points which passes through origin. 


calculated by taking into account the scattered neutron 
spectra and the number of neutrons of the input energy, 
convolved with the neutron detection response of the de¬ 
tector. 

The neutron transmission calculated above was used 
in the TOE data analysis in order to obtain the neutron 
spectrum at the BDS plane from the neutron spectrum at 
the scintillator plane, as shown in Fig. 3b. In this way, 
a direct comparison between neutron flux obtained by 
the scintillator and BDS can be acquired for an absolute 
calibration. A typical data set is shown in fig. 5a where 
the neutron spectra recorded from EJ232Q and EJ410 
scintillators are shown in comparison with the spectrum 
obtained by BDS, after applying a constant factor in or¬ 
der to match the neutron flux with that measured by the 
BDS. One can see a good agreement between the spectra 
in terms of spectral shape, which remains fairly consis¬ 
tent. 

The comparison between the neutron flux obtained 
from different scintillators with the one measured abso¬ 
lutely by the BDS over a number of shots is shown in 
figure 5b. In this case we have chosen the high energy 
range (2.5-20 MeV) for flux comparison in order to avoid 
discrepancies due to the down-scattered (lower energy) 
neutrons detected by the BDS, which is a time integrated 
spectrometer. As can be seen from the fig. 5b, the data 
points follow a linear trend over a range of neutron fluxes, 
the slope of which can be used as the final calibration fac¬ 
tor for the scintillators after applying other corrections, 
such as the scintillator detection response, the PMT gain 
and the transmission function for the setup. 


VI. CONCLUSION 

Employing a sub-ns neutron source driven by high 
power laser and absolutely calibrated, gamma insen¬ 
sitive bubble detector spectrometers, we have carried 


out the systematic calibration of three types of plas¬ 
tic scintillators namely EJ232Q, BC422Q and EJ410. 
Cross-calibration of the scintillators were done by imple¬ 
menting the neutron detection response (modeled semi- 
empirically for EJ232Q, BC422Q and using the one ob¬ 
tained from the manufacture for EJ410), PMT gain fac¬ 
tors obtained from the manufacturers and transmission 
of neutron through our experimental setup as simulated 
by FLUKA. Over a number of shots, spanning over a 
wide range of neutron fluxes between 10^ — 10^^ n/Sr for 
2.5-20 MeV neutrons, the comparison between the results 
obtained from the TOF data analysis and the BDS data 
follows a linear fit, which provides the final calibration 
factor for the scintillators. 
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